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Abstract 



■ The proton to proton polarization transfer coefficients K® , Kjj , and the 

proton to deuteron polarization transfer coefficients K®' , Kjj' , Kf, K^.' z ' , 
Ky Z ', K^ Z> , Ky' z ' and Ky X '~ yy have been measured in d(j?,J?)d and 



d(jf, ~ct)p reactions at E l £ h = 22.7 MeV, respectively. The data have been 



» ' z ' y y 
u(jf, ~ct)p reactions at E l ® 

compared to predictions of modern nuclear forces obtained by solving the 
three-nucleon Faddeev equations in momentum space. Realistic (semi) phe- 
"q ■ nomenological nucleon-nucleon potentials combined with model three-nucleon 

forces and modern chiral nuclear forces have been used. The AV18, CD Bonn, 
Nijm I and II nucleon-nucleon interactions have been applied alone or com- 
bined with the Tucson-Melbourne 99 three-nucleon force, adjusted separately 
^ ■ for each potential to reproduce the triton binding energy. For the AV18 poten- 

tial also the Urbana IX three-nucleon force have been used. In addition chiral 
NN potentials in the next-to-leading-order and chiral two- and three-nucleon 
forces in the next-to-next-to-leading-order have been applied. Only when 
three-nucleon forces are included a satisfactory description of all data results. 
For the chiral approach the restriction to the forces in the next-to-leading or- 
der is insufficient. Only when going over to the next-to-next-to- leading order 
one gets a satisfactory description of the data, similar to the one obtained 
with the (semi) phenomenological forces. 

21.45,+v, 24.70,+s, 25.10.+S, 25.40.Lw 
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I. INTRODUCTION 



A rich set of observables provided by three-nucleon (3N) reactions can be used to test 
modern nuclear forces. Presently there are two theoretical approaches to construct them. 
In the traditional approach the nucleon-nucleon (NN) potentials are derived in the frame- 
work of the meson-exchange picture alone or mixed with phenomenological assumptions. By 
adjusting parameters so called realistic, high precision interactions such as the phenomeno- 
logical AVI 8 [1] potential and the meson-theoretical CD Bonn [2] together with Nijm I and 
II [3] potentials are obtained. They provide a very good description of NN data below 
about 350 MeV nucleon laboratory energy. All these potentials have in common that they 
fit the large set of NN data with x 2 P er datum close to one, indicating essentially phase 
equivalence. 

In a more modern framework of chiral effective field theory, nuclear forces are linked 
to the underlying strong interaction between quarks and gluons. They are derived from 
the most general effective Lagrangian for pions and nucleons, which is consistent with the 
spontaneously broken approximate chiral symmetry of QCD, using chiral perturbation the- 
ory (xPT) [4]. The ;\PT approach gives a deeper understanding of nuclear forces than 
the traditional approach and allows to construct three- and more-nucleon forces consis- 
tent with the NN interactions [5-7]. In practice, various contributions to the nuclear force 
in the xPT framework are organized in terms of the expansion in Q/A, where Q is the 
soft scale corresponding to the nucleon external momenta and the pion mass and A is the 
hard scale associated with the chiral symmetry breaking scale or an ultraviolet cut-off. At 
present, the two-nucleon system has been studied up to next-to-next-to-next-to-leading 
order (NNNLO) in the chiral expansion [7,8]. At this order, the two-nucleon force receives 
the contributions from one-pion exchange, two-pion exchange at the two-loop level as well 
as three-pion exchange, which turn out to be numerically irrelevant. In addition, one has 
to take into account all possible short-range contact interactions with up to four derivatives 
and the appropriate isospin-breaking effects. In the three- and more-nucleon sectors, the 
calculations have so far only been performed up to next-to-next-to-leading order (NNLO) 
[7]. In this work, we will show the results corresponding to the latest version of the chiral 
NN forces introduced in references [9,8] and based on the spectral function regularization 
scheme [10]. 

Recent studies of few-nucleon bound states and of 3N reactions provided numerous in- 
dications that three-nucleon forces (3NFs) form an important component of the potential 
energy of three interacting nucleons [11-16]. In the traditional approach they are accounted 
for by adding model 3NFs, such as e.g. the 27r-exchange Tucson- Melbourne (TM) [17] 
or Urbana IX [18] interactions, with parameters adjusted to reproduce the experimental 
triton binding energy. Such a simple treatment allows to cure some of the discrepancies 
between data and theory [19-25]. In the approach based on x?T nonvanishing 3NFs appear 
in the next-to-next-to-leading order (NNLO) of the chiral expansion. In addition to the 
27r-exchange term two other topologies appear. One of them corresponds to a contact inter- 
action of three nucleons and the second to a contact interaction of two nucleons exchanging 
in addition one pion with the third nucleon. The two free parameters of these two terms 
are adjusted by fitting two independent 3N observables (e.g. the triton binding and the 
nd doublet scattering length). The quality of the description of the 3N observables is then 
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similar in both approaches [26]. The study of the details of the 3NF's is a lively topic of 
present day few-nucleon system studies. 

In the present paper we would like to analyze the proton to proton and the proton to 
deuteron spin transfer coefficients measured in d(~p*, J?)d and d(jf, d )p reactions, respec- 
tively, at E l p ab = 22.7 MeV [27-29]. The existing data for polarization transfer coefficients in 
elastic nucleon-deuteron (Nd) scattering are restricted to few experiments in the pd system 
( E iab = 1Q MeV [ 30 ] ; E l p ab = 19 MeV [31], E l p ab = 250 MeV [23], E l d ab = 270 MeV [32]) and 
to one measurement in the nd system (E l p ab = 19 MeV [33]). After a short description of the 
theoretical calculations in Section II in section III we show the data and compare them to 
different theoretical predictions. The summary and conclusions follow in Section IV. 

II. CALCULATIONS 

In this work we will employ two different methods to solve the nucleon-deuteron scatter- 
ing problem. Our first scheme is based on Faddeev equations. The nucleon-deuteron elastic 
scattering with neutron and protons interacting through a NN potential V and through a 
3NF V A is described in terms of a breakup operator T satisfying the Faddeev-type integral 
equation [34,15,35] 

T = tP + (1 + tG )\/ 4 (1) (l + P) + tPG T + (1 + tG )V^\l + P)G T. (1) 

The two-nucleon (2N) t-matrix t results from the interaction V through the Lippmann- 
Schwinger equation. The permutation operator P = P12-P23 + -P13-P23 is given in terms of 
the transposition P^ which interchanges nucleons i and j and Go is the free 3N propagator. 
Finally the operator v} 1 ^ appearing in Eq.(l) is part of the full 3NF V4 = vjp + v} 2 ^ + v} 3 ^ 
and is symmetric under exchange of nucleons 2 and 3. For instance, in the case of the ir — n 
exchange 3NF such a decomposition corresponds to the three possible choices of the nucleon 
which undergoes off-shell n — N scattering. It is understood that the operator T acts on 
the incoming state \<f> >= \qo > \<pd > which describes the free nucleon-deuteron motion 
with relative momentum q and the deuteron wave function \(pd >• The physical picture 
underlying Eq.(l) is revealed after iteration which leads to a multiple scattering series for 
T. 

The elastic Nd scattering transition operator U is given in terms of T by [34,15,35] 

U = PG 1 + PT + \/ 4 (1) (l + P) + V/ 4 (1) (l + P)G T. (2) 

We solve Eq.(l) in momentum space using a partial wave decomposition for each total 
angular momentum J and parity of the 3N system. To achieve converged results a sufficiently 
high number of partial waves have been used. Calculations with and without 3NF were 
performed including all 3N partial wave states with total two-body angular momenta up to 
j = 5. In the case when the 3NF is switched-off Eq.(l) is solved for J up to 25/2. When the 
shorter ranged 3NF is also active it is sufficient to go up to J < 13/2 only. In all calculations 
we neglect the total isospin T = 3/2 contribution in the 1 S'o state and use in this state a 
np form of the NN interaction. Such a restriction to the np force for the 1 Sq state does not 
have a significant effect on the polarization transfer coefficients. 
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The second scheme is based on the Kohn Variational principle, the ^-matrix elements 
corresponding to a 3N scattering state with total angular momentum J can be obtained as 
the stationary point of the functional 

[ J ^1 = J s s L i: + i{n SJ \H - E\n 8 -j) ■ (3) 

The wave function ^^sj describes a 3N scattering state in which asymptotically an ingoing 
nucleon is approaching the deuteron in a relative angular momentum L and total spin S. 
The parity of the state is given by (— 1) L . The wave function is expanded, using a partial 
wave decomposition, in terms of the pair correlated hyperspherical harmonic (PHH) basis 
as described in Ref. [25]. As in the Faddeev scheme, states up to J = 25/2 have been 
considered. 

In our Faddeev calculations the Coulomb interaction between two protons is totally 
neglected. A measurement of the neutron to neutron polarization transfer coefficients K% 
in neutron-deuteron (nd) elastic scattering [33] and their comparison to the corresponding 
pd data [31] shows that effects caused by the Coulomb force for this coefficient are non- 
negligible. These effects have been studied on a few polarization transfer coefficients [36] 
as well as in other polarization observables [25] using the Kohn variational principle in 
conjunction with the PHH basis. In these calculations the Coulomb force between the 
two protons has been considered without approximations and the results confirm sizable 
Coulomb-force effects in the energy range considered here. Therefore to remove possible 
Coulomb-force effects for the studied polarization transfer coefficients we proceed in the 
following manner. Using the PHH expansion we evaluate the studied polarization transfer 
coefficients without and with Coulomb force and employ the AV18 NN interaction. This 
is displayed in Figs. 1-3. In this manner we read off the shifts caused by the pp Coulomb 
force. Then we generate "nd" data by applying those shifts to our pd data. For the studied 
polarization transfers the Coulomb force effects are restricted mostly to forward angles and to 
the region around 9 cm pa 120°. For the proton to proton spin transfer coefficient K% they are 
of minor importance whereas for K% and K* the Coulomb force effects change significantly 
the magnitude of these coefficients around 6 cm pa 120° (see Fig. 1). For the proton to 
vector-deuteron spin transfers Coulomb force effects are rather small with exception of very 
forward angles, where they decrease significantly the magnitude (see Fig. 2). For K%' also 
some effects are seen around 9 cm pa 120°. In case of the proton to tensor-deuteron transfers 
shown in Fig. 3 only K% 2 (in the steep slope), K* z and K\ z exhibit large Coulomb force 
effects in the region of cm angles around 9 cm pa 120°. In the following figures we include 
both, the pd and "nd" data. 

III. RESULTS 

In Figs. 4-6 we show our data and compare them to theoretical predictions based on 
(semi)phenomenological NN potentials alone or combined with the TM99 [37] or Urbana 
IX [18] 3NFs, which have been obtained in the Faddeev approach, where we are able to em- 
ploy also non-local interactions. The corresponding comparison for chiral forces is presented 
in Figs. 7-9. Comparison with experiment always means the Coulomb corrected "nd" data. 

For the traditional approach based on high quality (semi)phenomenological interactions, 
we have taken the NN potentials AV18, CDBonn, Nijm I and II and combined each of them 



4 



with the 27r-exchange TM99 3NF, adjusting the cut-off parameter of TM99 individually 
to get the experimental triton binding energy. The resulting cut-offs for these potentials 
are respectively 5.215, 4.856, 5.120 and 5.072 (in units of the pion mass m n ). From the 
predictions of these potentials alone or combined with the TM99 3NF two bands, light and 
dark, respectively, were formed and shown in Figs. 4-6. 

For the proton to proton spin transfer coefficients (see Fig. 4) the effects of the Coulomb 
interaction are located in the region of cm. angles around 6 cm ~ 120°. In that region the 
differential cross section has its minimum. They are significant only for K y and K* and 
are practically negligible for K%' . The realistic potentials alone provide a good description 
only of and fail to reproduce the data for K y> and K*' , especially around 9 cm ks 120°. 
Including the TM99 3NF, and in case of AV18 also Urbana IX, changes only slightly the 
predictons for K% ■ For K y and K% the effects of these 3NFs are significant in the region of 
cm. angles around 6 cm m 120° and their inclusion leads to a good description of the data. 

For the corresponding proton to deuteron spin transfer coefficients (K% , K y and K* - 
see Fig. 5) effects of the Coulomb force are practically negligible at angles where data exist 
and they are seen only for . For these observables also effects of the TM99 and Urbana IX 
3NFs are small and the realistic NN potentials alone or combined with these 3NFs provide 
quite good description of the data. 

For the proton to tensor-deuteron spin transfer coefficients the Coulomb forces are sig- 
nificant for K y ' z ' ', Ky' z \ and K y ' z ' (see Fig. 6). In case of K z ' z ' and Ky' x '~ y ' y ' (see Fig. 6) 
they are negligible at angles where data exist. For these two coefficients also the effects 
of the TM99 and Urbana IX 3NFs are small and the NN potentials alone provide a good 
description of data. This is similar for K y z . For z and K y z the effects of these 3NFs 
are nonnegligible, especially in the region of angles around 9 cm ~ 120°. While for Ky' z ' the 
inclusion of the TM99 or Urbana IX 3NFs improves the description of the data, in case of 
K y ' z ' it shifts the theory away from the data points. 

Based on the chiral interactions we show in Figs. 7-9 two bands of predictions based on 
forces derived in NLO and NNLO. Each band is based on five predictions obtained with 
different cut-off combinations: (450,500), (600,500), (550,600), (450,700), and (600,700) (see 
[7] for more details). The results for the chiral NN potentials in NLO are shown by the light 
band. In NNLO the first time nonzero contributions from chiral three-nucleon interactions 
arise and the predictions based on the full chiral Hamiltonian in NNLO are shown by the 
dark band. 

It is clearly seen that the restriction to NLO only is quite insufficient even at low energy of 
our experiment. The predictions based on the chiral NN potential obtained in this low order 
are far away from the data. At NNLO one could show the NN force predictions alone. We 
refrain from doing that since this is ambiguous. As is well known unitarily transformed 2N 
forces, which do not affect two nucleon observables, lead to different results in the 3N system 
[38] . It is only the complete 3N Hamiltonian which provides unambiguous results for the 3N 
observables. Since the chiral approach systematically improves the nuclear force description 
with increasing order and provides strong internal links between 2N forces and forces beyond, 
we deviate here from the usual presentation exhibiting 3N force effects separately. This was 
done above in the standard approach since anyhow there is no internal consistency between 
NN and 3N forces. Now a view on Figs. 7-9 shows that the full NNLO predictions lead to 
a quite good description as the traditional approach including 3N forces. The exception are 
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for the proton to proton transfer, where the chiral approach differs from the data and 
Kf z ' for the proton to deuteron transfer, where it leads to an agreement with the data. In 
the standard approach it is opposite. 

It will be of great interest to see the outcome for NNNLO, where the NN forces have 
already been worked out [8] . At that order there contribute a whole host of parameter free 
3NFs including first relativistic effects. 

IV. SUMMARY AND CONCLUSIONS 

We presented new data for spin transfer coefficients in elastic pd scattering, both for 
the proton to proton and for the proton to deuteron transfers. They have been measured 
using a polarized proton beam with energy E l p ab = 22.7 MeV. The data have been corrected 
for Coulomb force effects using our theoretical framework of a hyperspherical expansion. 
This leads to "nd" data to which we compare our theory. In some cases the Coulomb force 
effects are quite significant, especially at 6 cm ~ 120°, where the differential cross section 
has its minimum. The theoretical predictions are obtained by solving the 3N Faddeev 
equations with two different dynamical inputs. One is the standard approach of so called 
high precision NN forces supplemented by the TM99 and Urbana IX 3NF. The other is 
an effective field theory approach constrained by chiral symmetry. In the first case where 
the forces have mostly phenomenological character we show NN force prediction separately 
in addition to the results obtained by adding the 3NFs. The inclusion of 3NFs clearly 
improves the description and the comparison with the data is quite successful. In case of 
the approach based on the chiral nuclear forces NLO is quite insufficient, but at NNLO 
the combined dynamics of NN and 3N forces does essentially equally well as the standard 
approach. Exceptions are the spin transfer coefficients from proton to proton and K y z z 
from proton to deuteron. The first is well described in the standard approach but not the 
second one. Just the opposite is true in the chiral approach. 
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[deg] 

FIG. 1. The nucleon to nucleon spin transfer coefficients in Nd elastic scattering at 
E% b = 22.7 MeV. The crosses are the pel data from [27-29]. The dashed line is the result of 
the hyperspherical harmonic expansion method with AV18 potential. The solid line is the corre- 
sponding result when the pp Coulomb force is included. 
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FIG. 3. The nucleon to deuteron spin transfer coefficients in Nd elastic scattering at 



Efa b = 22.7 MeV. The description of symbols and lines is the same as in Fig.l. 
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FIG. 4. The nucleoli to nucleon spin transfer coefficients in Nd elastic scattering at 
E? ah = 22.7 MeV. The crosses are the pd data from [27-29] and x-es are the corresponding "nd 
data" obtained by subtracting the effects due to the pp Coulomb force (see text for explanation). 
The light band results from theoretical predictions obtained with the AV18, CD Bonn, Nijm I and 
II NN potentials. The dark band is obtained when these interactions are combined with the TM99 
3NF. The solid line is the prediction of the AV18 + Urbana IX 3NF. 
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FIG. 6. The nucleon to deuteron spin transfer coefficients in Nd elastic scattering at 
22.7 MeV. The description of symbols, bands and lines is the same as in Fig.4. 
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FIG. 7. The nucleoli to nucleon spin transfer coefficients in Nd elastic scattering at 
^lab = 22-7 MeV. The description of symbols is the same as in Fig. 4. The light band results 
from theoretical predictions obtained with the NLO chiral potential with different cut-off parame- 
ters. The dark band results when in NNLO the NN- and 3N-forces are included. 
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FIG. 9. The nucleon to deuteron spin transfer coefficients in Nd elastic scattering at 



^lab = 22.7 MeV. The description of symbols and bands is the same as in Fig. 7. 
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